Charge stripe order from antiphase spin spirals in the spin-Fermion model 
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We revisit the ground state of the spin-Fermion model within a semiclassical approximation. We 
demonstrate that antiphase spin spirals, or 7r-spirals, whose chirality alternates between consecutive 
rows (or columns) of local moments, have, for sufficiently high carrier concentration, lower energy 
than the traditional Shraiman and Siggia spirals. Furthermore, 7r-spirals give rise to modulated 
hopping, anisotropic ID transport, and charge density wave formation. Finally, we discuss the 
relevance of 7r-spirals to the physics of charge stripe formation in cuprates, such as I^-zSr^CuCU. 

PACS numbers: 74.25.Ha, 75.10.Jm, 74.72.Dn 



I, INTRODUCTION 

The determination of the ground state of a strongly 
correlated system where conduction electrons interact 
with local magnetic moments is an issue that is rele- 
vant to many different areas of condensed matter physics. 
The usual starting point for such study is the so called 
Kondo lattice model (KLM), where the local moments 
are assumed to form a regular lattice. The KLM has be- 
come the paradigm for the study of the physics of heavy 
Fermion (HF) compounds and has been extensively in- 
vestigated both theoretically and with numerical meth- 
ods (for recent references see, for example, Refs. Q and 
0). When generalized to include a Heisenberg superex- 
change interaction, J, that favours antiferromagnetic or- 
dering among nearest neighbour local moments, the gen- 
eralized KML, or yet, the spin-Fermion (SF) model, has 
been argued to capture the essential physics of the high 
temperature superconductors (HTSC)£ 

Within the framework of the HTSC, the conduction 
electrons in the SF model are provided by the holes doped 
into the p x and p y orbitals of the ions, while the lo- 
cal moments are provided by the incomplete 3d x 2_ y 2 shell 
of the Cu ++ ions. Although this microscopic picture is by 
now well established, after 20 years of intense theoretical 
and numerical research^ the ground state of the HTSC 
is not yet completely understood. It is believed that the 
true nature of the ground state is determined by the com- 
petition between the different energy scales of the model: 
the hopping integral, t, the Kondo coupling, Jk, and the 
antiferromagnetic superexchange, J, as well as by disor- 
der. A large number of possible ground states have been 
obtained in the literature, including phase separation^ 
charge stripes formation^ and spiral magnetic order^ 
among others. 

From the experimental point of view, inelastic neu- 
tron scattering experiments within the superconducting 
phase of La2-£ C Sr 3; Cu04, x > 5.5%, revealed that dy- 
namical incommensurate (IC) spin correlations coexist 
with superconductivity^. In addition, it has also been 
observed experimentally the existence of IC "charge" 
peaks accompanying the IC magnetic order, with twice 
the incommensurability of the magnetic oneiAfi These re- 
sults were immediately and consistently interpreted as 



evidence of the formation of charge stripes that act as 
antiphase domain walls and apparently rule out other 
possible ground states such as Shraiman-Siggia spirals, 7 
which, despite being able to explain the IC magnetic or- 
der, require a uniformly charged ground stated 

In this brief report we revisit the ground state of the SF 
model for cuprate superconductors within a semiclassi- 
cal approximation. We demonstrate that antiphase spin 
spirals, or 7r-spirals, whose chirality alternates between 
consecutive rows (or columns) of local moments, have, 
for sufficiently high carrier concentration, lower energy 
than the traditional Shraiman and Siggia spirals. Fur- 
thermore, 7r-spirals give rise not only to IC magnetism 
but also to modulated hopping, anisotropic ID trans- 
port, and charge density wave (CDW) formation. Thus, 
our theoretical results indicate, for the first time, the pos- 
sibility of the coexistence between spiral magnetic order 
and CDW formation, in agreement with recent quantum 
Monte Carlo simulations for the SF model.-™ Finally, we 
discuss the relevance of 7r-spirals to the physics of charge 
stripe formation in the low temperature tetragonal phase 
of cuprate superconductors, such as I^-^Sr^CuO,^ 



II. THE SPIN-FERMION MODEL 

The Hamiltonian for the spin-fermion model can be 
written as H = H t + Hk + Hj where 

(i,j),a ia 

Hk = ~Y Y Si ■ p U aa < 3 Pi >p 

Hj = jJ^Si-Si' (1) 

p\ a creates a hole at the site i with spin projection a, t 
is the nearest-neighbour hopping integral, Jk is an an- 
tiferromagnetic exchange interaction between the Cu ++ 
local moment and the spin of the doped O -- holes, J 
is the Heisenberg antiferromagnetic superexchange, and 
[i is the chemical potential. In the large Jif/i limit, the 
above Hamiltonian is known to reduce to the t — J model. 
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Since conduction electrons and local moments in the 
SF model are independent degrees of freedom we can 
treat the local moments classically without affecting the 
mobile FermionsJ^ This is the basis of the semiclassical 
approximation used here. In this case, the Hamiltonian 
becomes quadratic in the Fermion fields and can be di- 
agonalized. For simplicity we choose a spin-quantization 
basis such that Si = (0, 0, S) at every site, and the 
Kondo term reduces to a shift on the chemical poten- 
tial for the different spin components of the doped holes, 

H K = (JK/2)SY,i(Pi,tPi,1 -PidPhl)- 

As it was originally proposed long ago by Shraiman 
and Siggiap for t > J (for I^-zSr^CuC^ i/J « 3) the 
hopping of the doped holes is favoured by a noncollinear 
configuration for the local moments of the spiral type (see 
also Ref . . Because of our choice for the local spin- 
quantization basis, as the spins spiral so does the local 
frame. This will affect the mobile holes that transform 
under the following SU(2) transformation 
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Here C is a unitary vector pointing to an arbitrary direc- 
tion in the plane perpendicular to the local spin, and a 
are the Pauli matrices. Under such SU(2) transformation 
the Fermion fields are then transformed as 
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FIG. 1: Inphase (top) and antiphase (bottom) spiral con- 
figurations for the magnetic ground state of the spin-Fermion 
model. For inphase (1, 0) spirals all spins rotate clockwise, for 
example, in every row m of copper ions. For antiphase (1,0) 
spirals, on contrary, the sense of rotation alternates between 
two consecutive rows, m and m + 1 for example. 
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and we usually write 



C/t(, 



°i,l 



q • Xj 



(3) 
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in such a way that the classical AF Neel state is 
parametrized by q = (tt,w) and has classical energy 
Eaf — —4NJS 2 , where N is the total number of local 
moments. 



III. INPHASE OR 0-SPIRALS 

Inphase or 0-spirals are parametrized as having an IC 
wave vector given by, for example, 



(5) 



This corresponds to Shraiman-Siggia original solution! 
(see Fig. and in this case the SF model reduces to 

ft=-*E co S ( q - (x ' 2 ' X3) )(4.^. a + fe.c.) 

<i,j),a V ' 

j. • ( 1 ' ( x * - x j) \ f -iC t ■ iC t \ 

_ i^sm K — J -L\(ie 4 cT T c i , i -ie l4 cT T c ia ), 



{id) 



Hj = -2NJS 2 {l-cosq x 



(6) 



The spiraling of the local moments favours the hopping 
of the conduction electrons in the direction of the pitch of 
the spiral via the H t term, at the price of some magnetic 
energy loss (Hj in Eq. © is larger than —4NJS ). In 
momentum space the above Hamiltonian reads 



k 

where 
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0.x 



6(k) = -2t 



sin fe T sm h sin fc„ 

2 u 



(7) 



(8) 



The above Hamiltonian can be diagonalized in momen- 
tum space and we obtain the dispersions 



e ± (k)=^ (k)± 



(9) 



with energy minima are located near ko = (7r/2, 7r/2) and 
symmetry related points in the magnetic Brillouin zone. 

We now calculate the total energy of the 0-spiral 
ground state. It is clear that, although the system pays 
magnetic energy in order to stabilize the 0-spiral, the gain 
in kinetic energy 
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FIG. 2: (Color online) Modulated hopping amplitudes for the 
two degenerate antiphase or 7r-spiral configurations at filling 
5 = 1/8. The blue (solid) curve is the hopping amplitude 
when the spins in the top row rotate clockwise (case drawn) 
while the red (dashed) curve is the hopping amplitude when 
the spins in the top row rotate counter-clockwise (not shown). 



is always larger than the magnetic energy loss, already for 
infinitesimal doping (as it has been demonstrated in Ref. 
\l^L a Dzyaloshinskii-Moriya term, not considered here, 
can shift this instability to finite doping, thus stabilizing 
the Neel ground state). Here we assumed that only the 
lower energy band in Eq. ^ is filled, due to the large 
gap (JkS)/2 (recall that Jx/t ^> 1), and we used that 
kp = V2ttS, where 5 gives the carrier concentration. 

The are two problems with the 0-spiral state: i) it 
does not break translational invariance in the charge sec- 
tor, as it is observed experimentally^ ii) it favours the 
transport in the direction parallel to the pitch of the 
spiral, while experiments have demonstrated that the 
anisotropic ID transport should be perpendicular to itil& 
As we shall demonstrate now, both issues are naturally 
incorporated by the antiphase or 7r-spiral state. 

IV. ANTIPHASE OR tt-SPIRALS 

Antiphase or 7r-spirals can be parametrized by the IC 
wave vector 

q w , = ((-) m flx,7r), (11) 



where m labels a certain row of Cu ++ ions, see Fig. 
The spin-Fermion model in this case is rewritten as (m 
and n label, respectively, rows and columns of local mo- 
ments) 



H t = -t cos Y( c (m,n+l) a C (r»,n)a + ^c) 

rn , n , a 

- t s[n ( n 1^( C l m +l,n) a C (mM) a + h.C.) 

m,n,a 

+ f^cos(ng 2; )(ie" I ' c c[ m+1 n)T c (min)1 + h.c.) 

m,n 

+ t^2cos(nq x )(ie 1(: cl m+ln)l c <m ^ + h.c.) 
_ JrS \r^ { t _ f \ 

Hj = 2NJS 2 cosq x . (12) 



There are important differences with respect to the 
case of the 0-spirals: i) in the 7r-spiral ground state, the 
spin independent part of the hopping perpendicular to 
the spiral pitch becomes modulated, see Fig. |21 ii) it be- 
comes maximal for parallel configurations of the local 
moments, as in positions n = 2,6; iii) it vanishes for 
antiparallel configurations for the local moments, as in 
positions n = 0,4,8. Furthermore, the modulation in 
the vertical hopping leads to a modulation on the carrier 
density, and, as a consequence, to a CDW instability of 
the stripe sort£ Finally, it is clear that the vertical hop- 
ping can be much larger than the horizontal one, thus 
exposing the ID nature of the transport in the 7t-spiral 
phase. 

We still have to show that the 7r-spiral phase has lower 
energy than the 0-spiral one. In momentum space the 
Hamiltonian reads 
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The above Hamiltonian can be diagonalized and the new dispersions are 



^(k) = ±jF I (k) ± 2/F>7(k), (14) 
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FIG. 3: Doping dependence of the energy of both inphase or 
0-spirals and antiphase or 7r-spirals. We see that although 
En > Eo at low doping, for larger 8 it becomes the true 
semiclassical magnetic ground state of the system. 



where 

Fj(k) = e 2 2 (k) + £ 2 (k) + 2A 2 (k) + {^f)\ (15) 
and 

FnQc) = (A 2 (k)-e 2 (k)e (k)) 2 +^) 2 (e 2 (k)+A 2 (k)). 

(16) 

Here 

A(k)=£sinfc y (17) 

is the CDW gap, with the property A(k + q) = — A(k), 
We see that Eqs. lfl4ll reduce to Eqs. ^ in the limit 



A = 0. For A ^ 0, in turn, both e\ +) (k) and e\ j (k) 
are always smaller than £q (k) , and thus 

E lin = iE{^ +) (k)eo,-4- +) (k)) 

k 

+ ^ £ (-)(k)e(/i- e (-)(k))|, (18) 

provides us with a rather large gain in kinetic energy. 
In particular, we found that for 5 « 5%, 7r-spirals al- 
ready have lower energy than 0-spirals, see Fig. We 
used units such that J — 1, with t — 3 and Jk = 5, 
which have the correct hierarchy tipically observed in su- 
perconducting cuprates (Jk > t > J). As a result, the 
stabilization of a striped CDW phase coexisting with ir- 
spirals breaks the translational symmetry in the charge 
sector, gives rise to IC magnetic correlations, and favours 
ID trasport perpendicular to the spiral pitch. 

V. CONCLUSIONS 

Using a semiclassical approximation we have revisited 
the ground state of the SF model. We have shown that 
7r-spirals have lower energy than 0-spirals and give rise 
to modulated hopping, anisotropic transport, and CDW 
formation. Thus, although such semiclassical analysis 
is, strictly speaking, only valid for large S, the possi- 
bility of the coexistence of spiral magnetic order and 
charge modulation is an appealing feature of the new 
semiclassical ground state here presented, captures the 
essential physics of the charge and magnetic IC order in 
La2- a; Sr 2 ;Cu04, and has been recently obtained numeri- 
cally with quantum Monte Carlo^ The stability of the 
semiclassical 7r-spiral ground state towards fluctuations 
has still to be demonstrated, but we believe that, as it 
happens for the case of 0-spirals in the t — J model^ 
the 7r-spiral phase can also be made robust through the 
inclusion of next-to-nearest neighbour hopping terms. 

The author acknowledges discussions with R. Doretto, 
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